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Abstract. We present the algorithmic details of the dynamical cluster approxima- 
tion (DCA) algorithm. The DCA is a fully-causal approach which systematically 
restores non-local correlations to the dynamical mean field approximation (DMFA). 
The DCA is in the thermodynamic limit and becomes exact for an infinite clus- 
ter size, while reducing to the DMFA for a cluster size of unity. Using the one- 
dimensional Hubbard Model as a non-trivial test of the method, we systematically 
compare the results of a quantum Monte Carlo (QMC) based DCA with those ob- 
tained by finite-size QMC simulations (FSS). We find that the single-particle Green 
function and the self-energy of the DCA and FSS approach the same limit as the 
system size is increased, but from complimentary directions. The utility of the DCA 
in addressing problems that have not been resolved by FSS is demonstrated. 



1 Introduction 

One of the most active subfields in condensed matter theory is the devel- 
opment of new algorithms to simulate the many-body problem. This inter- 
est is motivated by various physical phenomena, including high temperature 
superconductivity, magnetism, heavy fermions and the rich phenomenology 
occurring in quasi-one dimensional compounds. In the last few years, impor- 
tant progress has been made. Well-controlled results have been obtained by 
exact diagonalization and quantum Monte Carlo methods (QMC)|^. How- 
ever, these algorithms suffer from a common limitation in that the number 
of degrees of freedom grows rapidly with the lattice size. As a consequence, 
the calculations are restricted to relatively small systems. In most cases, the 
limited size of the system prohibits the study of the low-energy physics of 
these models. 

Recently, another route to quantum simulations has been proposed. Fol- 
lowing Metzner and Vollhardt[^ and Miiller-llartmann[0 who showed that in 
the limit of infinite dimensions, the many-body problem becomes purely lo- 
cal, a mapping to a self-consistent Anderson impurity problem was performed 
[ p^ . The availability of many techniques to solve the Anderson impurity 
Hamiltonian has led to a dramatic burst activity. However, when applied to 
systems in two or three dimensions this self-consistent approximation, re- 
ferred to as the dynamical mean field approximation (DMFA) , displays some 
limitations. Due to its local nature, the DMFA neglects spatial fluctuations 
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which are essential when the order parameter is non-local. This is the case in 
the cuprate or heavy fermion superconductors. Methods that at least include 
short-range spatial fluctuations are currently the object of intensive research. 
These techniques map the lattice problem to a self-consistently embedded 
finite-sized cluster, instead of a single impurity as in DMFA. 

The most promising of these techniques is the dynamical cluster approx- 
imation (DCA) The key idea of the DCA is to use the irreducible 
quantities (self energy, irreducible vertices) of the embedded cluster as an ap- 
proximation for the lattice quantities. These irreducible quantities are then 
applied to construct the lattice reducible quantities such as the Green func- 
tion or susceptibilities in the different channels. The cluster problem is solved 
by using a generalization of the Hirsch-Fye quantum Monte Carlo (QMC) 
method. 

In this paper, we discuss the novelty and the new perspectives offered 
by the DCA. We use the one-dimensional Hubbard model as a test case. 
We analyze the convergence of the DCA in comparison to that of the usual 
finite systems simulations (FSS). The choice of a one-dimensional model is 
particularly instructive. A success of the DCA in ID, where quantum fluctu- 
ations are the strongest, means that the method can work in any dimension. 
We show that the DCA produces results which are complimentary to those 
obtained by FSS. 




Fig. 1. Nc — 4 cluster cells (shown by different fill patterns) that partition the first 
Brillouin Zone (dashed line) . Each cell is centered on a cluster momentum K (filled 
circles). To construct the DCA cluster, we map a generic momentum in the zone 
such as k to the nearest cluster point K = M(k) so that k = k — K remains in the 
cell around K. 
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2 The Dynamical Cluster Approximation 
2.1 Formalism 

The DCA is based on the assumption that the lattice self energy is weakly 
momentuni dependent. This is equivalent to assuming that the dynamical 
intersite correlations have a short spatial range < L/2 where L is the linear 
dimension of the cluster. Then, according to Nyquist's sampling theorem]^ , 
to reproduce these correlations in the self energy, we only need to sample the 
reciprocal space at intervals of Ak ~ 27r/L. Therefore, we could approximate 
G(K + k) by G(K) within the cell of size (tt/L)^ (see. Fig. ||) centered 
on the cluster momentum K (wherever feasible, we suppress the frequency 
labels) and use this Green function to calculate the self energy. Knowledge 
of these Green functions on a finer scale in momentum is unnecessary, and 
may be discarded to reduce the complexity of the problem. Thus the cluster 
self energy can be constructed from the coarse-grained average of the single- 
particle Green function within the cell centered on the cluster momenta: 

^(K)^^EG(K+k) , (1) 

ic 

where N is the number of points of the lattice, Nc is the number of cells in 
the cluster, and the k summation runs over the momenta of the cell about the 
cluster momentum K (see. Fig. |l|). For short distances r < L/2 the Fourier 
transform of the Green function G{r) « G(r) + 0{{rAk)'^), so that short 
ranged correlations are reflected in the irreducible quantities constructed from 
G; whereas, longer ranged correlations r > L/2 are cut off by the finite size 
of the cluster. [|j 

This coarse graining procedure and the relationship of the DCA to the 
DMFA is illustrated by a microscopic diagrammatic derivation of the DCA. 
For Hubbard-like models, the properties of the bare vertex are completely 
characterized by the Laue function A which expresses the momentum con- 
servation at each vertex. In a conventional diagrammatic approach 

Z\(ki,k2,k3,k4) = ^exp [ir ■ (ki - k2 + ks - k4)] = iV^ki+k2,k3+k4 (2) 




Fig. 2. The DCA choice of the Laue function Eq. |j leads to the replacement of the 
lattice propagators G(ki), G(k2),--- by coarse grained propagators G(K), G(K'), 
... (Eq. y) in the internal legs of 'Pdca, illustrated for a second order diagram. 
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where ki and k2 (ka and k4) are the momenta entering (leaving) each vertex 
through its legs of G. However as £> ^ oo Miiller-Hartmann showed that the 
Lauc function reduces toQ] 

Zil3^oo(kl,k2,k3,k4) = l + . (3) 

The DMFA assumes the same Lauc function, AoMFAi}^i,^2,^7i,^4) — 1, 
even in the context of finite dimensions. Thus, the conservation of momentum 
at internal vertices is neglected. Therefore we may freely sum over the internal 
momenta at each vertex in the generating functional 'Pdmfa- This leads to a 
collapse of the momentum dependent contributions to the functional ^dmfa 
and only local terms remain. 

The DCA systematically restores the momentum conservation at internal 
vertices. The Brillouin-zone is divided into Nc — cells of size . 
Each cell is represented by a cluster momentum K in the center of the cell. We 
require that momentum conservation is (partially) observed for momentum 
transfers between cells, i.e., for momentum transfers larger than Ak ~ 27r/i, 
but neglected for momentmn transfers within a cell, i.e., less than Ak. This 
requirement can be established by using the Lauc function |^ 

^_DCA(ki, k2, ka, k4) = A^c5M(ki)+M(k3),M(k2)+M(k4) : (4) 

where M(k) is a function which maps k onto the momentum label K of 
the cell containing k (see. Fig. With this choice of the Laue function the 
momenta of each internal leg may be freely summed over the cell. This is 
illustrated for the second-order term in the generating functional in Fig. ^j. 
Thus, each internal leg G(ki) in a diagram is replaced by the coarse-grained 
Green function G(M(ki)), defined by Eq. 0. The diagrammatic sequences 
for the generating functional and its derivatives are unchanged; however, the 
complexity of the problem is greatly reduced since <^ N. We showed 
previously that the DCA estimate of the lattice free-energy is minimized 
by the approximation S(k) « i:(M(k)) + 0{Ak^), where S<Pdca/SG = S. 

One advantage of the DCA over other attempts to build self-consistent 
cluster techniques is that the DCA is fully causal. The spectral weight is con- 
served and the imaginary part of the single-particle retarded Green function 
and self-energy are negative definite. Hettler et al. used a geometric argument 
to derive a rigorous proof of the causality of the DCA [0 . 



2.2 One and two-particle quantities 

In the DMFA, after convergence, the local Green function of the lattice is 
identical to that of the impurity model. Though in the DCA, the coarse- 
grained Green function G{K) is equal to the cluster Green function, this 
quantity is not used as an approximation to the true lattice Green func- 
tion G{K). The correct procedure to calculate the lattice physical quantities 
within the DCA is to approximate the lattice irreducible quantities with those 
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of the cluster. The lattice reducible quantities are then deduced from the ir- 
reducible. In order to completely understand the DCA formalism, one must 
understand why reducible and irreducible quantities are treated differently. 
Consider a screened particle (a quasiparticle) propagating through the sys- 
tem, as pictured in Fig. ^. Here, we will assume that the screening length, 
e.g., due to Thomas-Fermi screening, ttf is short. This screening cloud is 
described by the single-particle self energy S{k^ iS) which itself may be con- 
sidered a functional of the interaction strength U and the single-particle 
propagator G(k, w). The different screening processes are described pertur- 
batively by a sum of self energy diagrams. Due to the short screening length, 
the propagators which describe these processes need only be accurate for 
lengths < Ttf- From the Fourier uncertainty principle, we know that the 
propagators at short lengths may be accurately described by a coarse sam- 
pling of the reciprocal space, with sampling rate Ak = -k /ttf- On the other 
hand, the phase accumulated as the particle propagates through the system 
is described by the Fourier transform of the single-particle Green function. 
Since this accumulated phase is crucial in the description of the quantum 
dynamics it is important that G(r) remain accurate at long length scales, so 
it should not be coarse-grained as described above. However it may be con- 
structed from the approximate self energy. Hence, the lattice Green function 
is given by, 

G(k) = . ^ v(M(v\\ ■ 

A similar procedure is used to construct the two-particle quantities needed 
to determine the phase diagram or the nature of the dominant fluctuations 
that can eventually destroy the single-particle quasi-particle. They are com- 
puted in the same fashion as the single-particle ones. The irreducible vertex 
r is the analogue of the self energy for two-particle quantities. In a general 
many-body theory, F has four entries describing the states of a pair before 
and after the interaction, 

r = r{kuk2,h,h) . (6) 




Fig. 3. Motion of a particle with its screening cloud 
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The diagonalization of F in spin space reduces to the usual density and mag- 
netic fluctuations in the particle-hole channel or singlet and triplet fluctuation 
in the particle-particle channel. For the sake of simplicity, we will use a generic 
notation below which describes any one of these four fluctuations. Momen- 
tum conservation ensures that there are only three independent momenta, 
r = r{q, k, k'). We assume, as for the self energy, that the irreducible vertex 
r = ^ is weakly momentum dependent. Hence, this quantity is equated 
to the corresponding cluster quantity within a cell surrounding the cluster 
momentum for each momentum involved in r{q, fc, k'). This gives, 

rQ+q(K + k,K' + k') =^q(K,K') . (7) 

The calculation procedure of the reducible two-particle quantities from the 
irreducible vertex F will now be described. The QMC cluster susceptibility 

Xc(g,i^,A") = (4+Q(T)cK'-Q(r')4(0W'(0)) , (8) 

calculated in the appropriate channel is the analogue of the single-particle 
Green function Gc- But this quantity is different from the coarse grained 
susceptibility x, 

XQ+q(K,K') = ^X;XQ+q(K + k,K' + k') , (9) 

k,k' 

where x is the lattice susceptibility. In the DCA, x and Xc are both related 
to r through the reduced Bethe-Salpeter equation, 

Xq+q = Xq+q + Xq+q^QXq+q , (10) 

represents either Xc or x- Xq+q is obtained from the single-particle Green 
function as follows, 

X^+q = '5K,K'^^G(K + k)G(K + k + Q + q) . (11) 

In Eq. ^ above, the sum is restricted to cluster momenta. This constitutes 
a significant simplification of the lattice problem which may be otherwise 
intractable. Yet, as the cluster size increases, the full solution of this reduced 
Bethe-Salpeter equation may require a significant amount of computer stor- 
age and CPU time. We thus restrict ourselves to momenta Q where a given 
type of fluctuation is likely to diverge, i.e., Q = for the density, singlet or 
Q = (tt, tt) magnetic fluctuations. The inversion of each of the two equations 
then yields 

= Xc' - x'.^' + ■ (12) 

The charge (ch) and spin (sp) susceptibilities \ch,sp 

are deduced from x 
Xch.sp{c^,it^,i) ^ ^^^^2 X! ACTCT/Xq(K, icJm;K',iw,„/) , (13) 

^ KK'mm'crcr' 

where \aa' — 1 for the charge channel and Ao-cr' = crcr' foi" the spin channel. 
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2.3 The DCA algorithm 



Table 1. Steps of the DCA algorithm 

0. set r = 

1. calculate the coarse-grained Green function G 

2. calculate = + E 

3. using Hirsch-Fye calculate the cluster Green function Gc 

4. compute a new estimate of 17 = — G^^ 

5. repeat 1-4 until S converges 

6. accumulate bins of measurements on Gc, Xc 

7. calculate the lattice Green function G, spectral weight, susceptibilities... 



The DCA iteration procedure is set forth in Table ^ It is started by setting 
the initial self energy to zero. This self energy is then used to compute the 
coarse-grained Green function G'(K). The latter is used to compute the host 
Green function Q{K.)^^ = G(K)~^ + ^(K) which must be introduced to 
avoid over-counting diagrams. Q serves as the input to the QMC simulation 
to yield a new estimate of the cluster self energy. The procedure is repeated 
until S converges. This typically happens in less than ten iterations. The 
number of iterations decreases when Nc increases since the coupling to the 
host is smaller for larger clusters {F sa ©(l/iVc)) H). The convergence test is 
made on the ratio p, 

\T.n^old{i^n)\ 

Once convergence is reached, the remaining single and two-particle measure- 
ments are made in a final QMC iteration. As in a usual QMC simulation, 
bins of measurements are accumulated as discussed in the next section. One 
should, however, bear in mind that these measurements are performed with 
G instead of the true lattice Green function G. Hence, the determination of 
the lattice quantities requires an additional step in which the coarse-graining 
equation is inverted. This is done in a separate program. 



3 The Hirsch-Fye QMC algorithm 

The cluster problem is solved using one of the available finite system tech- 
niques. The Hirsch-Fye (HF) algorithm which was originally proposed 
for impurity problems, is easily adapted to embedded cluster simulations as 
required in the DCA. We recall here the main steps of its derivation. For 
details we refer the reader to the original Hirsch and Fye paper or the review 
by Georges et al [T^ . One starts by writing the partition function in the 
Grand-canonical ensemble. 
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Then the exponential is broken up into Nr imaginary-time sUces and the 
Trotter approximation is employed to separate the kinetic, K, and interacting, 
V, parts of the Hamiltonian H, 



Nr Nr 
1=1 1=1 

where Nt-A^ — (3. The error made during the last step is of the order of 
(Z\t)^. The next step is to replace this interacting problem by one consisting 
of non- interacting particles moving in a fluctuating field Xa , where the index 
a stands for space and imaginary-time coordinates. This is done by employing 
the Hirsch-Hubbard-Stratonovitch (HHS) transformation, 

p[-zl,t/(n„T-i)(««i-|)] _ ^ " ^[Aa;„(n„T-«ci)] 



Xa,=±l 

After the HHS transformation, the quartic operator V becomes quadratic 
like K. In this new form a trace over the fermion degrees of freedom can be 
exactly performed. It yields, 

^ = E n det[{G:{x^))-^] , (18) 

where the sum is over the configurations of HHS fields Xa- The dimension 
of the Green function matrix G% is [Nt x Nc)"^. It is at this point that 
the HF algorithm differs from the Blanckenbecler-Sugar-Scalapino (BSS) 
algorithm In the HF, the simulation deals directly with the matrix G^. 
In the BSS the locality of the action in time is used to reduce the matrix-size 
in the simulation. However, this cannot be done when action is non-local in 
time, as in this case where the cluster is coupled to a host. 

The algorithm starts with an initial Green function Gc — G and a cor- 
responding initial configuration of the fields Xa- It then sweeps the space- 
time lattice by proposing fiips in the fields Xa x'^ — —Xa- The heat-bath 
or the Metropolis algorithm is used to determine if the change will be ac- 
cepted. The probability P{xa) of the configuration Xa is proportional to 
nCT=Ti det[{G'^{xa))~^]- The transition probability is 

_ Y{.=ud^A{G':{x^))-'] 

Y{.=ud^i[iGt{x^)y^] • ^ ^ 

If the new configuration is accepted, the Green function is updated by using 
the relation 

The initial field configuration is chosen with all x^ — 0. The above equa- 
tion is then used to construct a physically realistic field configuration (i.e.. 
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Xa = 1 for all fields). The system is then warmed up by sequentially stepping 
through the space-time lattice, proposing changes at each space-time site 
Xa — > —Xa- The change is accepted if the transition probability is greater 
than a random number between and 1. Typically, the warm up phase lasts 
for about a hundred space-time sweeps before measurements begin. It is nec- 
essary to perform a few complete space-time sweeps in order to produce more- 
or-less independent measurements. For clusters, the Hirsch-Fye algorithm is 
very efficient and stable at low temperatures. 



4 Numerical Aspects 

One difficulty encountered with the DCA algorithm is that a reliable trans- 
form from imaginary-time quantities, in the QMC part, to Matsubara fre- 
quencies, for the coarse-graining part is needed. A careful treatment of the 
frequency summation or the imaginary-time integration is crucial in order to 
ensure the accuracy and the stability of the algorithm and to maintain the 
correct high-frequency behavior of the Green functions. We need to evaluate 
the following integral 

f0 

Gc{K,iu;n)= / dre^'^"Ge(K,T) . (21) 
Jo 

But from the QMC, we know the function GcCK, t) only at a discrete subset of 
the interval [0,/?]. As it may be readily seen by discretizing the above equa- 
tion, the estimation of Gc(K,ia;„) becomes inaccurate at high-frequencies. 




Fig. 4. (a) Imaginary-time Green function at the Fermi point for U = 2, (3 = 5 for 
FSS (solid lines) for iVc=4,8,16 (from top to bottom) and for DCA (dotted lines) 

for A^c=4,8,12 (from bottom to top), (b) The corresponding spectral weights for 
FSS (solid lines) A^c=4,8,12 (increasing value at w = 0) and DCA (dotted lines) 
iVc=4,8,12 (decreasing value at a; = 0). 
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This is formalized by Nyquist's theorem which tells us that above the fre- 
quency Wc = unpredictable results are produced by conventional quadra- 
ture techniques. A straightforward way to cure this problem may be to in- 
crease the size of the set of r-points where the Green function is evaluated. 
But, this renders the QMC simulation rapidly intractable as seen in the pre- 
vious section. A much more economic way to avoid the problem is to use an 
approximate method that is asymptotically exact. 

Second-order perturbation theory is enough to obtain the correct asymp- 
totic behavior. For instance, we compute the Matsubara-frequency Green 
function from the imaginary-time QMC Green function as follows [Eol 



Ge(K,zc^„) = Gept(K,zc^„)+ / dre^"^ (G,(K, t) - G,p* (K, t)) 

Jo 



(22) 



The integral is computed by first splining the difference Gc(K, T)—Gcpti^; ''')> 
and then integrating the spline (a technique often called oversampling). 

Once convergence is reached, G — Gc, and the QMC Green function Gc 
may be analytically continued using the Maximum-entropy method (MEM). 
Unfortunately, there is no reliable way to perform the direct analytic contin- 
uation of i7(K). Pade approximants lead to very unstable spectra because 
of the QMC statistical noise contained in Z'(K). The binned imaginary-time 
Green function data accumulated from the cluster calculation must be used to 
obtain lattice spectra from which S(K) may be deduced. To obtain the self- 
energy and spectral- weight function A(k, lo) of the lattice in real frequencies, 
we first compute the cluster spectral-weight A('K,lu) by using the Maximum 





Fig. 5. (a) Imaginary-time Green function at the Fermi point for 17 = 4, /3 = 5 for 
FSS (solid lines) for 7Vc=4,8,16 (from top to bottom) and for DCA (dotted lines) 
for 7Vc=4,8,12 (from bottom to top), (b) The corresponding spectral weights for 
FSS (solid lines) A'^c=4,8,12 (increasing value at uj — 0) and DCA (dotted lines) 
A^c=4,8,12 (decreasing value at lj = 0). 
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entropy method |^ for the inversion of the foUowing integral equation 

G{K,t) = J dcuj^^AiK^u) , (23) 

where G'(K, r) is the imaginary-time Green function obtained from the QMC 
simulation of the cluster. From ^(KjLj), we obtain the coarse-grained Green 
function in real frequencies by Kramers-Kronig analysis. Finally, we solve 
equation (1) to obtain Z'(K,cj) from G(K,uj). 



5 Example: the One-dimensional Hubbard Model 

In this section, we apply the DCA to the one-dimensional Hubbard model at 
half-filling. Such a test is interesting for several reasons. First, this test will 
help verify whether recent applications of the DCA in two dimensions are 
reasonably anticipated to be accurate. In one dimension, quantum fluctua- 
tions are stronger than in higher dimensions. Hence, one intuitively expects 
the DCA to be less efficient in one than in higher dimensions. So, if the DCA 
accurately captures the physics in one dimension, then it is highly likely to 
capture the physics of two and three dimensions accurately. Second, although 
the DCA is known to become exact in the limit of an infinite cluster, an ex- 
tensive, systematic analysis of the convergence of the DCA has not yet been 




0.2 0.4 0.6 0.8 1 -5 -4 -3 -2 -1 1 2 3 4 5 
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Fig. 6. (a) Imaginary-time Green function at the Fermi point for U = 6, /3 = 5 for 
FSS (solid lines) for A'^c=4,8,16 (from top to bottom) and for DCA (dotted lines) for 
A^c=4,8,12 (from bottom to top), (b) The corresponding spectral weights for FSS 
(solid lines) Afc=4,8,12 (broad to narrow gap) and DCA (dotted lines) Nc^4,8,12 
(narrow to broad gap). Note that once a gap opens in the spectral function, the 
maximum-entropy analytic continuation procedure becomes unreliable. Q This un- 
reliability is believed to be the source of the non-systematic nature of the peaks in 
the spectral functions. Although unreliable for fine structure, the qualitative feature 
of the existence of a gap is accurately depicted in this figure. 
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performed. Finally, it is desirable to illustrate differences between the DCA 
and a widely-applied FSS method in a well-studied, non-trivial problem. 

The ID Hubbard Hamiltonian is widely accepted as the most relevant 
model in the high temperature regime of the ID organic materials, for which 
a significant volume of experimental data is available. It is well-studied and 
provides a non-trivial test of the DCA. The Hubbard Hamiltonian reads: 

H = -tJ2ictc^+la+hc) + UJ2r^^1n^l . (24) 

i i 

with a next-nearest-neighbor hopping and an on-site repulsion t and U , re- 
spectively, on a one-dimensional lattice. We set t = 1 throughout this study 
and measure all energies in terms of t. We work at half-filling, where the QMC 
is free of the fermion sign problem, eliminating one possible source of errors 
in both the FSS and the DCA. This allows us to easily isolate actual discrep- 
ancies between the FSS and the DCA. Although the ID Hubbard model can 
be solved exactly, the FSS QMC is currently the only reliable method used 
to compute finite temperature dynamics. 

We now turn to a comparison of the imaginary-time Green function at the 
Fermi-point of the DCA and of FSS. We see that the two methods converge 
differently to the exact result. The quantity 

G{k„m) = Ud.^^^ , (25) 

is useful to measure the strength of correlations. At the Fermi-point and at 
a fixed temperature T, it varies from — i for J7 = 0, an uncorrelated system, 
to for C/ = cxo, a highly-correlated system. In Fig. ^(a), ||(a), and ^(a), 
we display G(t) for C/ = 2, 4, 6; = 4, 8, 16; and /3 = 5 for FSS and for the 
DCA. These parameters were chosen to illustrate the three situations that are 




Fig. 7. Real (a) and Imaginary (b) part of E at the Fermi point for U = 2, P — 5 
for FSS (solid lines) for iVc=4,8,16 and for DCA (dotted lines). 
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generic in the model: the metalhc, the pseudogap, and the insulating regimes, 
respectively. In all cases, G{kp, (3/2) decreases with increasing cluster size for 
FSS while for the DCA it increases. 

This behavior marks the fundamental difference between FSS and the 
DCA. At low temperatures, in FSS, the correlation length is greater than 
the lattice size. Thus, the effects of correlations are overestimated for smaller 
clusters because the systems are artificially closer to criticality than a system 
in the thermodynamic limit. This tendency is reduced by increasing the clus- 
ter size, which moves the system in the direction of the thermodynamic limit. 
The situation is radically different in the DCA where the system is already in 
the thermodynamic limit. The DCA approximation restricts correlations to 
within the cluster length. As the cluster size increases, longer range correla- 
tions are progressively included. Thus, the effects of the correlations increase 
with the size of the cluster. 

This analysis is supported by the spectra shown in Fig. ||(b), ||(b), and 
^(b). For [/ = 2, FSS shows a pseudogap for Nc = 4. This pseudogap vanishes 
at Nc = 8, 16 to become a peak at w = 0. The peak corresponding to Nc = 16 
is sharper that of Nc = 8. In contrast, the DCA starts with a sharp peak 
at Nc ~ 4. This peak progressively broadens as Nc is increased to 8 and 16. 
The FSS and the DCA peaks seem to converge to the same limit consistent 
with the results of G(t). For [/ = 4, the system presents a pseudogap. The 
convergence to this pseudogap is in conformity with the above analysis. The 
FSS system goes from a gap to a pseudogap when Nc is increased from 4 
to 16. The DCA evolves from a central peak at A^c = 4 to a pseudogap at 
Nc — 16. When [/ = 6, the system is gapped; the FSS shows a large gap at 
Nc — ^- This gap decreases for Nc — ^ and Nc = 16. 

In Fig. ^(a), ||(a), and ^(a), we show the real part of the self energy 
&i (3 — h for various interaction strengths, U. For U — 2 (Fig. |^), where a 




Fig. 8. Real (a) and Imaginary (b) part of E at the Fermi point for [/ = 4, /3 = 5 
for FSS (solid lines) for Afc=4,8,16 and for DCA (dotted lines). 
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pseudogap exists for Nc = 4, Ui via FSS has two solutions at non-zero lu for 
the equation — — = 0. This happens when the slope '^^'^^^^ is 

greater than unity. This derivative decreases and becomes smaller than 1 for 
Nc = 8, 16 corresponding to the single peak in A{(jj). In the DC A, 
slowly increases with Nc but always remains smaller than 1. 

The Imaginary part of the self energy (w) shown in Fig. |^(b), ^(b), and 
^(b), also has a monotonic behavior in which the limiting value is bracketed 
by FSS and DC A results. A general trend that emerges is that S2Ujj) has a 
local minimum in the vicinity of a; = in the metallic regime (Fig. 0)j which 
minimum deepens sharply when a pseudogap appears (Fig. and finally 
E2(lo = 0) vanishes in the insulating regime (Fig. ^). 



6 Summary 

We have extensively analyzed the DC A by comparing it to FSS. A coherent 
picture emerges from the investigation of the single-particle properties of 
the ID Hubbard model. By systematically underestimating the effects of 
correlations, the DCA converges to the limiting value of a given quantity 
from a starting point that is opposite to the starting point of FSS. Thus, the 
DCA can help resolve situations where one is unable to draw a conclusion 
from FSS alone. For example, in spite of considerable work with FSS, it has 
not been possible to make a definitive statement as to the existence of a 
pseudogap in the weak-coupling (U smaller than the band-width) regime of 
the two-dimensional Hubbard model at half-filling. It was impossible to get 
a definite answer from FSS because, as expected, the pseudogap seemed to 
vanish when the cluster size was increased iQ. Recently, using DCA, Huscroft 
et al. showed that there is a pseudogap for 4 x 4 to 8 x 8 clusters |l^ . From the 




Fig. 9. Real (a) and Imaginary (b) part of E at the Fermi point for [/ = 6, /3 = 5 
for FSS (solid lines) for Afc=4,8,16 and for DCA (dotted line). 
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analysis above, since the DCA underestimates the effects of correlations, this 
demonstrates unambiguously the existence of the pseudogap in the half-filled 
Hubbard Model in the thermodynamic limit. 

The DCA has shown great promise and, even though a young method, 
has already illustrated its great utility. Still, much work remains to be done. 
One immediate application is the the nature of the ground state of the two di- 
mensional Hubbard model; is it a Mott insulator or antiferromagnetic (Slater) 
insulator? Another is the problem raised by Anderson concerning the low en- 
ergy behavior of 2D coupled doped Hubbard chains; are they Luttinger or 
Fermi liquids? Both of these problems have proven difficult to tackle 
using FSS approaches alone. Finally, we have not discussed the convergence 
of the two-particle quantities. A systematic comparison with FSS as in this 
work will be the subject of a future study. 
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